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Fig. 1. - Notation for the lectures. 



1. Quantum Information Processing with linear òptics 

Thcrc arc at lcast three ways to undcrstand these lectures on linear òptics quantum 
computing: 
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• these are engineering lectures on how to build an optical quantum computer. We 
will give a step by step approach on how to use beam splitters, phase shifters, singlc 
photon sources and detectors to efficicntly simulate an ideal quantum computer. 
In this sense, these lectures can really be thought of as being an engineering recipe, 

• they can be thought of as lectures on quantum error correction. In order to use the 
simple elements mentioned above to build quantum information processing devices 
we will have to rely heavily on ideas and concepts of quantum error correction. By 
putting elements together and using detectors, a simple error model will emerge 
(projection of a qubit in the Z basis) and we will use fault tolerant error correction 
to ensure that the quantum information is protected. We will even calculate a value 
for the aceuracy threshold, 

• or these can be seen as lectures on what quantum information is. These lecture will 
teach you that maybe we do not understand what quantum information really is. 
Some years ago it was thought that a quantum computation would be implcmented 
by doing very precise unitary transformations. Here we will do very non-unitary 
operations (projections) on physical qubits but we will still be able to show that at 
a higher level this can mimic a unitary evolution efhciently. Linear òptics will allow 
us to do a set of transformation that can be simulatcd efficicntly by a classical com- 
puter. It is only when we condition our operations on the results of measurements 
that the system becomes hard to be simulated classically. 

IT. Quantum òptics and quantum information. - These lectures will be divided into 
three sections. First an introduction to quantum òptics which will include some historical 
notes on the use of quantum òptics for quantum information processing (QIP). The 
second section will focus on using linear òptics with singlo photon sources and detectors 
(sometimes called the LOQC or KLM theory). The third section will focus on quantum 
error correction and LOQC. 

1"2. Quantum Computation. - Quantum information promises to solve some problems 
that have no known efficient classical algorithms. This is done by harnessing the size of 
Hilbert space and the ability to manipulate the quantum state. 

David DiVincenzo pQ has suggested a set of five criteria that are sufficient for a device 
to be used as a quantum computer: 

1. A scalablc physical system with well characterised qubits. 

2. The ability to initialise the state of the qubits to a simple fiducial state, such as 
|000...). 

3. A "universal" set of quantum gates such as genèric one qubit gates and a two qubit 
gate, 

4. A qubit-specific measurement capability. 

5. Long relevant decoherence times, much longer than the gate operation time. 
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l'3. Why òptics?. - Many physical systems have be proposed for quantum computa- 
tion. Each has its own advantages and disadvantages. For example, in an ion trap, two 
qubit gates are relatively casy to implcmcnt. Howevcr. isolating the ions from the envi- 
ronment is difficult, the ions' motion being vulnerable to decoherence. In NMR quantum 
computing, the most successful proposal to date, nuclear spin states can be handled di- 
rectly. However, the ability to initialise a simple fiducial state is difhcult. Also, many 
nuclear spins must be present in order for a sizable signal to be observed. 

The first proposal for quantum computing involved single photons as qubits [7] . One 
advantage for using single photons in QIP is that quantum òptics is a well developed ficld. 
Photons decohere slowly, an ideal trait for quantum computing. Photons also travel well, 
one reason why they are so widely used for communication. Another advantage is that 
photons can be experimentcd with at room temperaturc. 

One of the major obstacles for using single photons for QIP is that photons do not 
interact directly, making two qubit gates very difhcult. Mode matching with single pho- 
ton is also a problem. Another difhculty with using single photons is how to produce and 
detect them. To date the best single photon sources and detectors both have cfhcicncics 
well below any known thrcshold for sealable quantum computation. 

l'4. Quantum Òptics. - 

l'4.1. Classical Electromagnètic Field. When we deal with electromagnètic waves 
classically wc look for solutions to the source free Maxwell eqs. 

(1) V 2 E(r,í)-^E(r,í) = 
The solution have a plus and minus frequeney part [21 E| : 

(2) E(r,í) = i ^(^) 2 [a fc u,(r)e-^ í + <.u*(r)e^ t ] 

fe ^ ' 

The Ufe(r) are orthogonal mode functions, usually plane wave mode functions ( Ufc(r) = 
^7= exp(zk • r). where è' A ) is the unit polarisation vector and V is the volume), and the 
a,k arc dimensionless (complex) amplitudes. 

The energy of a classical electromagnètic field is given by 



(3) H= I^(E 2 +B 2 )dr = £(^ja fc < 

l'4.2. Quantisc. Wc quantise the electromagnètic field by turning the coefficient dk 
into operators and imposing the commutation relations: 



(4) [à i ,à]] = s iJ 

H>àj] = [à i , uyj 



(5) [à i ,à j } = [àt,à\} = o. 
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In terms of these operators the energy is given by 
(6) H = ^2,hw k {à\à k + | 



k 



If we let 



the energy of the electromagnètic field becomes 

ü P 2 _l 1 2-2 
ii = — míü x 

2m 2 

from which we can recognise the harmònic oscillator. 

We call à an annihilation operator and à' a creation operator, the reason for this will 
be shown in eqs. ijSJl and ©. 

The eigenstates of H are labeled \n) and we call them Fock states. Wc define the 
number operator n by n\n) — à^à\n) — n\n). 

What is the cffcct of à on \n)l We can answer this by finding the number of photons 
in à\n): 



AtA2|„\ 



So à\n) is a Fock state with n — 1 photons. We define this to be ^4|n — 1). In a similar 
way we can show that à'\n) is a Fock state with n + 1 photons, _B|n + 1). But what is 
the form of A and B? 

We can use the fact that (n\h\n) = n: 

(n\n\rt) = (n\à'à\n) — A(n\a í \n — 1) 

= (n - í\A 2 \n - 1) = A 2 = n 
(n\n\n) = (n|à'ò|n) = (n\àà^ — l\n — 1) 

= B 2 - 1 = n 



So we see that à acts as an annihilation operator and à' as a creation operator: 

(8) à f |n) = Vn + ï\n + 1) 

(9) à\n) = - 1) 

from which we deduce that n is a semi-positivc number. 

For n — 0, a|0) = and H\0) = ^fuv\0} representing vacuum fluctuation energy of 
the lowest eigenstate of the Hamiltonian. 



Linear Òptics Quantum Computation: an Overview 



7 



The Fock states form an orthonormal set: 

(n\m) = S nm 

OO 



?7=0 



And each Fock state may be built up from creation operators: 



(10) 



Exercise: show that In) 



|0) 



|0). The matrix form of the annihilation operator 



can be seen by using the identity operator twice: 

OC OO 

à = Ml = |n)(ri|à 



Tl = 



m— 



E VnTT\n)(n + 1| 



n=0 



In matrix form: 



/ s/ï ••• \ 
V2 ••• 
\/3 ••• 



V 



/ 



In eq. Q) we wrote à in terms of x and p. It is often more convenient to write the 
annihilation operator as a linear combination of the two Hermitian operators: 

Q1+1Q2 



where Qi (equivalent to x from before) corresponds to the in-phase component of the 
elèctric field amplitude of the spatial-temporal mode and Q2 (equivalent to p from before) 
corresponds to the out-of-phase component. 

This gives the commutation relation [Qi^Q/l = 2i and the Heiscnberg uncertainty 
relation 

AQ1AQ2 > 1 



With the annihilation and creation operators defined in this way it is instructive to 
look at the position and momentum wave functions for Fock states. First note that 
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Í'(Qi) = (Qi\i>) and ^(Qa) = (QalV')) where we know that Qi|Qi) = Qi|Qi) and 
Q2IQ2) = Q2IQ2) ■ Also note the usual normalisation rules: {Q\\Qx) = 6(Qi — and 
(Q2IQ2) — S(Q2 — Q 2 )i ^(^ — x ) being thc Dirac Delta function. 

Whcn we look at the position and momentum wave functions for the vacuum statc 
we have 0] 



(Qi\o)=MQi) = 

(Q 2 |0) = MQ 2 ) = 




And in general for a Fock state \n) we have 



(Qi\n) = Vn(Ql) = 




whcre H n (x) is the Hermite polynomials. 

l'4.3. Minimum Uncertainty States. Thc Hciscnberg Uncertainty Principle: for two 
non-commuting observables A and B, the uncertainty in each is given by the relation: 

AÀAB>±\<[À,è]>\ 

wherc (Ai) 2 = (i 2 ) - (i) 2 . 

For example, consider the case when A = x and B = p. The uncertainty principle 
becomes the well known result: AxAp > | 

A coherent state |a) is a minimum uncertainty state. We can define it as the eigen- 
state of an annihilation operator: à\a) = a\a) Alternatively, we can define it with the 
displacement operator D(a): 

\a) = D(a)\0) = cxp(aàï - a*à)\0) 
= exp(-^)exp(aàt)|0) 

Exercise: show the last line for the above coherent state definition is correct using 
the Campbcll-Bakcr-Hausdorff operator identity: 

e A+È = e A e È e -í [A,B] 

for the case when [À, [À, B]] = [È, [À, È]] = 0. 
Exercise: show that àD(a)\0) = cü£'(q;)|0). 

If we work out AQ1AQ2 for a coherent state we find that AQi = AQ 2 = 1, showing 
that coherent states are minimum uncertainty states, since AQ1AQ2 = 1- 
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Coherent states are not orthogonal since 

((3\a) =exp{~(\a\ 2 + \f3\ 2 )+af3*} 
and as such form an over complete set 

J \a){a\d 2 a = 7r 

If we look at the wave functions for a coherent state ((Qi\a) = ïp a (Qi) and (Q2\a) = 
ï>a(Q2)) |H we find that 

tp a (Qi) = exp 
ï ) a(Q2) = tt~* exp 

where we have decomposed the complex amplitude a into its real and imaginary parts: 
a = \ {qo + ipo)- From this we see that a coherent state is a Gaussian in Qi, Q2 phase 
space. 

Squeezcd states are another type of minimum uncertainty state. Squeezed state have 
one of AQi or AQ2 "squeezed" below 1 at the expense of the other: either AQi < 1 < 
AQ2 or AQ2 < 1 < AQi such that AQ1AQ2 — 1 is always true. 

We can definc squeezed states with the Squcczc operator 

5( e ) = cx P (i e *à 2 -i e (at) 2 ) 

where e — re ltp . 

If we define a rotated complex plane: P1+ÍP2 = (Qi+iQ2)e~^' p then the uncertainties 
can be written as APi = e~ r and AP 2 = e r . Squeezed coherent states are defined to be 
Ke)=£>(a)S(e)|0>. 

If we look at the wave functions for a squeezed coherent state 0] wesee that (Qi\D(a)S(e)\0): 

_i § / 2e(Qi - lo) 2 , «Vi n i ^ 
ir 4e2exp(^-e ze + -Q2Q1 - gPo<7oJ 

where we have again decomposed the complex amplitude a into its real and imaginary 
parts: a — 5(90 + *Po) ■ As with the coherent state, we see that a squeezed coherent state 
is a Gaussian in Qi, Q2 phase space. 



(Qi ~ Qo) , i 



P0Q1 - zPoQo 



(Q2 - Po) 



- 4<7o<?2 + gPoqo 
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l'5. Linear Òptics. - What does "Linear Òptics" mean? 

An optical component is said to be linear if its output modes b. are a linear combi- 
nation of its input modes aj: 

k 

Linear optical components are madc up of phasc sl·iiftcrs and bcam splittcrs. A phasc 
shifter is defined by the transformation U(P ( f > ) : \n) — > e m ^\n). That is, (ój)"|0) — > 
e m *(àl) n \0). 

The Hamiltonian for a phasc shifter is given by the number operator: U = e l ^ H = 
e uj>n _ e ï0a l a_ Tj S mg this see wc that e l ^° ta |n) = e m ^\n) Exercise: show that the phase 
operator takes a coherent state \a) to |e^a). 

A beam splitter is defined by the transformation matrix 



(11) U 



-e'^sini 
cos 9 



where the input modes are related to the output modes via aj\0) — > X)m^mi a m|0) 
Exercise: The Hamiltonian for a beam splitter can be given by H = e^à^b + e~ l ^àw , 
check that this is Hcrmitian. 

The unitary transformation for a beam splitter then looks likc expí— i6{e í ^ct) k ài + 

e àkà\)^j , where k and l are the two modes being acted upon by the beam splitter. 

Consider we have the state \mn)\2 incident on a beam splitter. The output would 
look likc: 

à{cos6» + 4e" í0 sin6» ( a}e" í0 sin 9 + à\ cos 6) |00) 



vrrúrú 



For examplc: 



|10) -> cos6»|10) +e-^sin(9|01) 
|01) -> -e^sinfl|10) + cos0|OT) 
|11) -> -\/2e # cos0sin0|2O) + (cos 2 9 - sin 2 0)|11) 

+y/2e-^ cos é> sin 6» 1 02) 
|20) -> cos 2 (9|20) + V2e~^cos0sin0|ll) + e" 210 sin 2 6>|02) 
|02) e 2 ^sin 2 0|2O) - V2e^ cos 6» sin 0| 11) + cos 2 9\02) 



Exercise: confirm these examples. 
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l'6. Previous Suggestions With Òptics. - Photons are an ideal means for communi- 
cation, they travel fast and only interact weakly with their surrounding and thus have a 
long dccoherence time. The problcm with using solcly òptics for quantum computation 
is how to get photons to interact. It is relatively easy to achieve single qubit gates, but 
how do we establish two qubit gates? This isn't the only problem, we also have to achieve 
two qubit gates that are cfRcient in resources so our system will be scalable. One way to 
make photons interact is to induce a cross phase modulation between two optical modes. 
One example of this is the Kerr effect. 

l'6.1. Quantum Optical Fredkin Gate. The first proposal for quantum computation 
was the quantum optical Fredkin gate [ül {7\ . The gate was realised with single photon 
òptics using the Kerr Effect. A Fredkin gate is a three qubit gate that acts as a controlled 
swap. Its truth table is given below. 



Cj o>i c Qjq b 





1 1 

10 10 

110 11 

10 10 

10 1110 
110 10 1 

111111 



Table I. - The lògic for a Fredkin gate. Here subscript i refers to the input qubits, o refers to 
the output qubits and c is the control qubit. 

In this scheme, logical corresponds to the vacuum mode |0) and logical 1 is the 1 
photon Fock statc 

The optical Kerr effect is defined by a material with an intensity dependent refractive 
index, that is, a nonlincar crystal that has an index of refraction n proportional to the 
total intensity /: n — n + n 2 E 2 — n + n 2 I, where n is the normal refractive index and 
n.2 is the correction term necessary for Kerr materials. 

The Kerr Hamiltonian is of the form 

Hi = —hxà\àia2à2 

where \ i s a coupling constant which depends upon the third-order non-linear suscepti- 
bility for the optical Kerr effect and à\ and à 2 are the annihilation modes for the input 
light. 

The optical setup for the Fredkin gate first proposed in [7j is shown below in fig. [21 
In this setup we have a Mach-Zchnder interferometer with a Kerr media in cither arm. 
In the top arm the Kerr media has the control beam incident on it. 

The beam splitter on the left side, BS1, has 8 = j, <fi — and the beam splitter on the 
right hand side, BS2, has 8 = f , <j> = tt- The effect of the Kerr media is to act the unitary 
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operator U = exp {—iíHj\ = exp (ieà 1i àc í c) in the top arm and U = exp (ieb^bd^ ct) in 
the bottom arm, whcrc e = t\. 

For example, consider the input state \cab, d) = |010,0) 

|010, 0) ->-^(|010,0) + |001,0)) 

^ e itóteat V e£tyt<! (|010,0) + [001, 0» 
v2 

= -^( e ieeteata |010,0) +e" 8tytd "|001,0)) 
V 2 

= 4=(|oio,o) + |001,0>) -> |010,0) 
v 2 

where mode c? is always in the vacuum mode |0) and remember that 
Ncxt consider the input |101,0) 

|101,0) - -J= (| 101, 0)-| 110,0)) 



-^(-e íeeteats |110,0) +e íeStSd " td *| 101,0)) 
v2 

-L(-e"|110,0) + |101,0)) 

i(-e le |110,0) + e íe |101,0) + |101,0) + |110,0». 
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If we choose e = ir this gives 1 110, 0). With this choice of e we can show that all the lògic 
gates are as in the Fredkin gate lògic table. 

There are two major problems with this scheme: (1) it is difficult to achieve the 
high non-lincarities, cspccially thosc rcquired for e = 7r phase change; (2) at high non- 
linearities the crystal exhibits other dctrimental affccts, such as absorption. 

16.2. Cavity Quantum Electrodynamics. Using non-linear crystals is not the only 
method to induce a coupling between optical modes. In [H] Turchette et. al. used cavity 
quantum electrodynamics to induce a two qubit interaction between photons. To do this 
single àtoms in a Fabry-Perot cavity were used. 

We start with the Jaynes-Cummings Hamiltonian for the interaction of photons with 
a two level atom: 

(12) H = —hüJo(-Z) + hutà^à + g(í$(j- + àcr + ) 

where lüq is the frequeney of the atom, lo the frequeney of the field and g is the coupling 
constant between the field and atom. 

The first term in the above Hamiltonian is the atom energy term iï a tom = — \hu)çjZ. 
We have a two level atom and we say the energy of the top level (|1)) is and 
that of the bottom level |0) is — ^hu>o. From the diffcrence in energy we have H at0 m = 
ÍHu) (\ï)(l\ - |0)(0|) = given that Z = -(|1)(1| - |0)(0|). The # field = fiiüà^à 

term comes from the harmònic oscillator describing light (and we have dropped the \ fno 
term as it plays no role here) Finally, there is the interaction term Hi which shows that 
an atom in the state |1) can decay to give a photon, or a photon can excite the atom from 
a state |0). a± are the usual raising and lowering operators, a+ = |1)(0| = h(X + iY\, 
(X- = |0)(1| = \{X — iY). The terms àV + or àcr_ arc dropped from this Hamiltonian 
by the so called rotating wave approximation. 

We can rewrite the Jaynes-Cummings Hamiltonian in the following form, as done 

(13) H = tkvN + 5Z + g(iïa- + àcr+) 

where N = à^à + \Z and 5 = \h(üj$ — un . We may neglcct from any further analysis 
as it leads to a constant phase factor that we can neglect. We thus have a Hamiltonian 
that is reduced to: 

(14) H = 5Z + g(àï<7_ +àcr+) 

The goal is to utilise a single atom to obtain a non-linear interaction between two 
photons. The quantum information will be represented by orthogonal polarisation photon 
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states, using a cavity with a single atom to provide the non-linear interaction between 
the photons. We can modify the Hamiltonian (|14(l to incorporate two photons of slightly 
difïercnt freqüències and take photon polarisation into account. A thrcc lcvel atom that 
is excited to a diffcrent cncrgy lcvel (of approximately the same energy) by the two 
orthogonal polarisation, a so called V type atom, sufnces for this. 

Turchcttc at. al. uscd Cesium àtoms and two photons, one circular ly polariscd (mode 
à) and the other of slightly diffcrent frequeney also circularly polarised (mode ò), both 
generated from weak làser pulses. The lògic was represented by the two orthogonal 
polarisation modes |l ± ) a |l ± )b. The phases induced by interaction with the atom were 
then measured using hetcrodyne measurement. 

The following logical transformations were measured: 



|l-) |l->6 "> |l-)a|l->6 

|l + )a|l-)6-e^|l + )a|l-) fa 

in a |i + ) 6 - e *ii-) Q |i + ) b 

|1 + >«|1% ^ e < <*·+** +A >|l + ) |l + )i 



where 4> a ps (17.5 ± 1)°, <p h ps (12.5 ± 1)° and A ps (16 ± 3)°. 

The bcauty of this interaction is that |1 + 1 + ) gocs to e (' Pa+I ' ,t+A ) |1+1+), wherc intu- 

itively one might expect that |1 + 1 + ) would go to e( iPa+Vb ^ . This is the photon 

non-linearity. If A were 0, the system would be linear and the photons would not inter- 
act. This is the key for implementing quantum lògic with thesc mcthods. Unfortunatcly 
it is mudi harder to get higher non linearities such as A = 90° and does not appears to 
be a promising approach when compounding many gates. 

l'7. Progress with Linear Òptics. - Another direction to devise quantum information 
processor is to focus on linear òptics alone. In this section we show what progress was 
made in this direction before the proposal by Knill et. al. UJ. 

l'7.1. Decomposition of unitàries. Reck et. al. showed that we can break any unitary 
into a set of linear optical components |lüj . The unitary transformation considered is 
acting on the creation operators, so one is not able to immcdiately apply this method to 
say to the CNOT unitary. 

A unitary operator U of dimension N can be decomposed as follows: 

(15) U = ÍTjv,iV-l ' TV,iV-2 ' ïjv,JV-3 ' • ' ■ Tjf-l,N-2 ■ TV-l,iV-3 • • • ÏV-DJ 

where T p q is the N dimensional identity with the {p, q} elements replaced with the beam 
splitter matrix l|ll|) and D is an N x N matrix with phases on the diagonal. 

The general linear optical network for a unitary matrix U is a triangular array of 
beam splitters and phase shifters, shown in fig. [3J 



15 



Mirror 




Fig. 3. - General linear optical network for a unitary matrix U. 



As an cxample, what does the linear optical circuit for the following unitary look likc? 




Exercise: check that this is unitary. 

Here we have a unitary matrix with N=3, so we can say that 

u= (t 3i2 -T3, 1 ·t 2)1 · j d)~ 1 = 

We find that 

fi 

T 3 , 2 = COS 01 

V -e* sinó'l 



^ J 2,l J 3,l J 3,; 





1 sin 6»i 



COS#! 



Ï3.1 = 



COS 6*2 



-e* siné» 2 




1 




g-#2 g in( 


COS 6*2 



2,1 



COS #3 

-e* sin 6» 3 




s"^ 3 sin 6» 3 



COS 6*3 






1 
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and 



D = 








Exercise: show that a solution to this is given by the angles: 9\ = 12.8°, </>i = n, 
9 2 = 20.4°, 2 = 0, 9 3 = 63.8, <fe = 0, </> 4 = rr, 5 = 0, 6 = 0. 

Exercise: What transformation does this implement on the input state (a|0) +/3|1) + 
7|2))|10)? What is the transformation if we measure 1 10} in modes 2 and 3? This will 
be important latcr! 

l'7.2. Optical Simulation of Quantum Logic. In Cerf et. al. proposed a scheme 
for quantum lògic with only linear optical devices and a single photon. To simulate n 
qubits a single photon is put into 2™ different paths. 



A VNOT gate is given by a beam splitter with 9 = ? and 



VNOf =±( 1 1 



whcrc 



V2 V 



|0) 9 = |01)^-^(z|10) + |01)) 
|l) 3 = |10)-i=(|10)+i|01}) 

Thrce simple gate implcmentations proposed by Cerf et. al. are given in hg. 

W ,/ v (6) 



|1) 



|1> 



10 > |o> 



|o> 



11) 



11) 



T ^ 

^ 

m i- 



Fig. 4. - A Hadamard, a CNOT and a reverse CNOT. 



With these a universal set of gates is possible. 
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To implement a Hadamard gate we use a 8 = ^ , <j) = — -| beam splitter and two — ^ 
phasc shiftcrs, as in part (a) of fig. 01 

|0) fl = |01)-»-i|01)-^~(|01)+i|10)) 

_^(|oi> + iH)|io)) = ^(|oi) + |io)) 

|l) g = |10)^|10)^i=(i|01} + |10}) 

-,-L(i|oi) + -»|io)) = -J=(|oi)-|io)) 

To implement a CNOT we encode a qubit in position and polarisation. The location 
is the control and the polarisation the target. For the control: \0) q = 1 01) , |1) 9 = 1 10) . 
For the target: |0) g = \H), \l) q = \V). The circuit for this is a polarisation rotator on the 
upper arm, as in part (b) of fig. If a photon is present in the top arm, its polarisation 
will be flipped: 

( a|01}+/3|10» <g> (j\H) + 6\V)) 

-> a-/\01)\H) + a5\01)\V) + f3j\lO)\V) + f3S\10)\H) 

To implement a reverse CNOT we simply need a polarising beam splitter (PBS), 
where horizontal is reflected. As before, the location is the control and the polarisation 
the target. This is in part (c) of fig. 0] 

( a|01)+/?|10» ® (j\H)+S\V)) 
-> a7|01)|ír) +aó\10)\V) +/3j\10)\H) + /3S\01)\V) 

The problem with this scheme is that n qubits requires 2 n paths which in turn re- 
quires 2" — 1 beam splitters to setup. This is not scalable. With one qubit encoded in 
polarisation we still need 2"" 1 optical paths. 



2. Linear Òptics Quantum Compution 

In this section we will dcvelop ideas and concepts for linear òptics quantum computing. 
We will show that using beam splitters, phase shifters, single photon sources and detectors 
we will be able to simulate efficiently, i.e. using a polynomial amount of resources, an ideal 
quantum computer. It is known that a device using only linear òptics (beam splitters, 
phase shifters) is not more powerful than a classical computer. We will see that once 
we use the result of measurement and feed it back to future linear gates, we can now 
simulate a quantum computer. 
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2T. Assumptions in LOQC. - Before we start explaining the details of linear òptics 
quantum computing (LOQC), we first recali what "linear" means when we say linear 
òptics. 

For an optical element to be linear, each output mode &t is a linear sum of the input 
modes à\: 



where M is a unitary matrix. 

We use beam splitters and phase shifter to implcment linear òptics: 

A phase shifter is described via: U = e m *, where n is the number of photons in the 

mode. 







Fig. 5. - A phase shifter of phase 



A beam splitter is described by: 



(16) 



U 



cos 



sm( 



cos( 



where the input modes arc related to the output modes via 



4lo)^^[w4jo) 



e * 



Fig. 6. - A beam splitter of angle 9 and phase 



The bàsic resources that are necessary for linear optical quantum computing arc singlc 
photon sources and detectors that can distinguish between 0,1 and 2 photons. 
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2'2. Qubits in LOQC. - In LOQC we encode qubits with dual rail lògic. That is, lògic 
that is encoded using the physical location of a photon. Logical is given by |0) q = |01) a 6 
and logical 1 by \l} q = |10) o 6, whcre mode a is one spatial mode and mode b anothcr. 

We could just as easily use photon polarisation: |0) g = \H), \l) q — \V). We can 
change the representation of the qubit from polarisation to spatial mode using a polarising 
rotator and a polarising beam splitter as can be seen bclow in hg. [7] 



a\H)+P\V) 



Q|01> ai , + /3|10) ai) 




a\H) + p\V) 



Fig. 7. - Transforming from dual rail to polarisation encoding. 



The polarisation encoding of qubits is useful to understand some of the experimental 
implcmcntation of LOQC but for this and subsequent lectures we will be using dual rail 
encoding for qubits. 

2'3. Qubit Operations . - For a physical system to be a viable candidatc for quantum 
computation, we need a universal set of gates, as stated in point 3 of DiVincenzo's five 
criteria pQ, seen in section ll'21 Such a universal set is comprised of arbitrary single qubit 
gates and the CNOT operation |13| . 

Wc will show that single qubit operations are easily performed using linear elements 
and almost always error free. However, as we saw last lecture, it is quite difHcult to 
make photons interact. We will show that with the use of single photon ancillas and 
photo-detection we can make a two qubit gate (a control phase gate). The non-linearity 
will still be there, in the form of measurement. Unfortunately this two qubit gate will 
not be unitary but using quantum teleportation and quantum error correction we will be 
able to approximatc a unitary operation using only a polynomial amount of resources. 

2'4. Single qubit gates. - All single qubit gates can be implcmented with just beam 
splitters and phase shifters. To see this, note that any single qubit unitary, U, can be 
decomposed into rotations about the Z and Y axis in the Bloch sphere as follows 

U = é a R z (P)R y ( 1 )R z (5) 

where R z {6) = e'^"* and R y (<p) = e"^. Exercise: show that R x (d) = R z {^)R y {-d)R z { 
If we can show how to make arbitrary rotations about the Z and Y axis we can 

perform any arbitrary single qubit operation. 

Rotations about the Z axis can simply be performed using a phase shifter on the top 

mode of a qubit: 
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1 

l'O 



o — 



Fig. 



Performs a rotation of (j> about the Z axis: Rz (</>). 



Considcr the cfFect of the above circuit on the qubit a|0) g + 



a\0) q + p\l) q = ot|01> +j9|10> 

-► a|01) +/3e^|10) 
= e^ /2 (e-^ /2 a|0) 9 + e^ /2 /3|l) 9 ) 
= e^ 2 e-^/ 2 (a|0) 9 +/3|l) ? ) 
= e^ 2 iïz(0)(«|O) 9 +/?|l) 9 ) 



And we see that, up to an irrelevant global phase, a Rz(4>) has been performed. 

Rotations of —26 about the Y axis require a beam splitter of angle 8 and angle = 0: 



e d> 



Fig. 9. - Performs a rotation of —20 about the Y axis: Ry(— 2$). 



We can see this by first remembering the effect of a beam splitter on a photon in 
modes 1 and 2: 



— ► (cos(6*)q| + sin^àj) 
à\ — > (— sin(#)àj + cos^àíj) 
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We can consider the effect of the above circuit on the qubit a|0) q + 



a\0) q + f3\l) q = a\01) +/3|10) 

-> a(cos(0)|Ol) - sin(0)|lO}) + /3(cos(0)|lO) + sm(0)|Ol}) 
= cos(0)(a|Ol) + /3]10>) - sin(0)(a|lO) - /3|01>) 



«(alO^+fll),) 
ifr(-20)(a|O),+j9|l),) 



2'5. Two guòií gates. - Now we have to address the question of how to make a CNOT 
gate with just single photon ancillas and photo-detection. We will actually be looking 
at how we can perform a CSign, the relation between the two is apparent from the fact 
that HZH = X. The CSign gate has the action \x) L \y) L -► e™ x -»\x) h \y) L , x,y € {0, 1}. 



H 



X 



H 



Fig. 10. - The equivalence between a CNOT and a CSign. 



We will first show how to perform a probabilistic CSign gate on photonic qubits. 
From here we will show how the CSign probability of success can be boosted arbitrarily 
close to 1 using teleportation. In the next section we will see how this could be simplificd 
even further by using quantum error correction. 

To perform a CSign gate wc nccd a bàsic transformation: ct|0) + /3\1) + 7|2) — > 
a\0) + /3|1) — j\2). This is called a non-linear sign shift NS-i. Exercise: why is the 
transformation non-linear? 

2'5.1. Nonlinear sign shift gate. As the above transformation is non-linear we cannot 
use linear òptics alone. We find that with the use of ancilla modes and photo-detection 
we can indeed perform this transformation using two extra modes, one with a singlo 
photon and the other in vacuum as shown in the circuit below. The input state looks 
like (a|0>! + 7 |2)i) ® |l)a ® |0) 3 . 
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l¥> 



Fig. 11. - The NS-! gate. 



Where 1 = 22.5°, 0i = 0°, 2 = 65.5302°, 2 = 0°, 6 3 = -22.5°, 3 = 0° and 

<f> 4 = 7T°. 

The action of this circuit is to apply the following unitary matrix on the input modes: 




Does this look familiar? Remember lecture 1. 

Whcn we measure a single photon in mode 2 and vacuum in mode 3 we have the statc 
a|0)i + j8|l)i — t 1 2) i output in mode 1. The probability of measuring a single photon in 
mode 2 and vacuum in mode 3 is ^. 

Exercise: show that this circuit indeed performs the transformation ct|0) + /3|1) + 
-f\2) — ► a\0) + (3\1) —7|2) when mode 2 and 3 are measured in |10) and that this succeeds 
with probability 1/4. 

2'5.2. Controlled Sign gate. With the use of two NS-i gates we can make a CSign 
gate, as seen in the fig. below. 

Lets work through a general cxamplc, say we start with the input statc Q\ = a\0) q + 
/?|l)q = a|01) +/3|10) in modes 1 and 2 and Q 2 = -f\0) q + 5\l) q = 7 |01) + S\10) in modes 
3 and 4. 
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NS.-, 
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i 


/ 


45 




-45° 



Fig. 12. - A probabilistic CSign gate. The probability of success 



Aftcr thc first bcam splittcr wc havc 



a 7 |0101) + a<5|0110) + /3 7 |1001) + /3<5|1010) 
«710101} - ^11100) + ^10110) + ^11001) 

+ ^|0011) -^|2000) + ^|0020) 



Qi 



Q 2 







NS.., 


3 


i 


/ 


45 


3 


-45° 



Fig. 13. - Step 1 of a probabilistic CSign gate. 
After the two NS-i gates we have 



oí8 oió 3^ 

«7|0101> - t=|1100) + ^=|0110> + ^=11001) 

+^|0011) + ^|2000) -^|0020) 
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Qi 



Q 2 



V 




NS.-, 


3 




/ 


45 


3 


-45° 



Fig. 14. - Step 2 of a probabilistic CSign gatc. 



After the final bcam splitter we have 



a7|0101) + ad\01W) + /3 7 |1001) - /3«5|1010) 
a 7 |00) ? + oí|01), + fa\W) q - (3S\ll) q 



Qi 



Q 2 







NS.., 




À 


/ 


45 





-45° 



Fig. 15. - Step 3 of a probabilistic CSign gate. 

This gate works with a probability of j^. That is, there is a probability of jg that 
we will measure the ancilla state 1 10) in both NS-i gates. It is important to notice that 
when the result of the measurement of the ancilla state is other than |10), the gate fails 
but we know that this failure has occurred. 

It is also interesting that during the gate, between the ±45° beam splitters, the 
photons are not in the dual rail lògic encoding but rccombine at the last beam splitter 
to get back into the right encoding. Thus we go out of the qubits Hilbert space to make 
this gate and go back at its output. 

The state 



(17) 



i ( |0101) + |0110) + |1001) - |1010) 
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is the first building block to performing a CSign on photonic qubits with probability 
arbitrarily close to 1. 

2'5.3. Tclcporting qubits through a gatc. What wc have shown so far is how to 
perform a probabilistic CSign gate on photonic qubits. This is not cnough to allow 
scalable quantum computation. The probabilistic CSign gate will be used as an entangled 
state production stage. 

To perform a two qubit gate we could teleport ^3] the qubits through the gate as 
first shown by Gottesman and Chuang |15j . 

Say we have two arbitrary qubits \ipi)q and \tp2)q that we want to perform a CSign 
on. We could teleport both thesc qubits and then apply a CSign on them, as scen in the 
fig. below: 




Fig. 16. - An application a CSign between qubits |V>i)g and \tp2) q - The dashed-dotted box is 
the state preparation area. In this case just the Bell states |£?oo) need to be prepared. 

Exercise: given that \B m ) q = ^(|00) + \i>i) q = c \0) q + Cl \l) q , \ip 2 ) q = 

do\0) q + di|l) q (all logical qubits, |co| 2 + |ci| 2 = |cío| 2 + \d\\ 2 = 1), and that we measure 
in the logical basis, show that this circuit gives a CSign between \ipi) q and \ip2)q- 

Since CSign applied twice gives the identity, we can add two CSign's to fig. as 
seen in fig. El 
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^1)9 { 




Fig. 17. - An application a CSign between qubits \ij)i)q and \tp2)q 



We now have a CSign inside the dashed-dotted box, the state preparation area. Now 
we want to look at the effect of conjugating the "corrections" with the CSign's. To do 
this we will nccd the idcntity: 







/ 1 

10 
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V -1 / 
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CSign ( a x ® 1 ) CSign = a x g) a z 
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^ 




í° 





1 


\ 





1 
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-1 
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X 



Fig. 18. - Shows the identity CSign ( a x ® 1 ) CSign = a x <g) a z 
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and thc identities 



CSign^ll g) cr x J CSign = a z ® a x 
CSign (ct 2 <g> ij CSign = a z ® 1 
CSign (]1 ®a z ) CSign = 1 ® <r 2 



Applying these identities gives the following circuit: 

IV'l), { = 



IV'2), < 




Fig. 19. - An application a CSign between qubits \ipi) q and \ip2)^ 



The problem of performing a CSign has now been turned into a state preparation 
problcm. We can apply the CSign gate (fig. I12fl probabilistically offline. Once we know 
we have a successful CSign state (eq. i)17f) '). we can proceed and teleport thc information 
qubits, as in fig. I19I This way we do not corrupt the quantum information. 

2'5.4. Teleporting with the C-Sign entangled states. If a gate fails with probability 
/ and we perform error correction, after we do k of these gates we will have a success 
probability of / . With teleportation, constructing two qubit gates reduces to a state 
preparation problem. Thc critical element is that we know when the state preparation 
fails and thus we can ensure that it does not corrupt the quantum information. To 
succeed we need to attempt the state preparation j times to implcmcnt k gates. 

The next stcp is to show how we perform telcportation with linear òptics? 



2'5.5. Basic Teleportation with linear òptics. It is known |12| that we cannot distin- 
guish bctwccn thc four Bell states with linear òptics alone. With a single beam splitter 
the best that we can do is distinguish bctwecn thc Bell states, and thus do teleportation, 
with a success probability of i. We will show how this is done and then generalise the 
telcportation so that the probability of success is as near to one as desired. Although 
we cannot achieve a probability of one as this would requirc perfect Bell state distin- 
guishability, wc can use quantum error correction and limit the left over error so that it 
is smaller than what is required from the aceuracy threshold theorem. 
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Consider teleporting the state a|01) + /5 1 10) . We need only look at teleporting the 
first mode of this state: 



I . 1 ) \ 



(Ri 



(R2)-i !?>„> 



c, 



Fig. 20. - Telelportation with linear òptics. 



First the entangled resource state ^(|01)34 + 110)34) is made, shown below. At this 
point the state is 



( a|01)i 2 + /3|10>i 3 )- ? =(|01>34 + |10> 34 



1 / 



= -^=(a|0101) 123 4 + a|0110)i 2 34 + /3|1001)i 2 34 + /3|1010)i234 



(R1H 



c, 



Fig. 21. - Telelportation with linear òptics: making the entanglement resource state. 
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After the second beam splitter we have: 

-> - (\/2a|0101) - a|1100) + a|0110) + /3|1001) + /3|0011) - /?|2000) + /3|0020)) 
Remembering that the second mode was unchanged here. 




Fig. 22. - Telelportation with linear òptics: performing a partial Bell measurement. 

If the detectors measure a total of 1 photon (R1+R2 = 1), we can recover the original 
information in modes 4 and 2: 

i?! = 0,i? 2 = l: |(ai|01>+j9|10» 
i? 1 = l,i?2 = 0: |(-a|01>+>9|10» 




Fig. 23. - Telelportation with linear òptics: photo-detection. 
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In the case that we measure R\ = 1, R2 = we need to perform a phase shift correction 
gate Ci and we have seen how to do that in the previous section. The probability of 
successful teleportation is \ + \ = 4. If the detectors measure or 2 photons the 
teleportation fails by collapsing the state in the mode |0) or |1) and this is equivalent to 
projecting the qubit in the Z basis, as will be explaincd in the next section. This occurs 
with a probability of \. 




Fig. 24. - Telelportation with linear òptics: applying a correction. 

2'5.6. The telcported C-Sign. Now we can piece togethcr what we have worked out 
so far to makc a CSign gate. First we have the CSign entangled state production step: 




Fig. 25. - CSign state production 
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Fig. 26. - Teleportation. 



The success probability for the state production stage (fig. is jg, however this is 
donc offiinc. This state preparation fails (on average 16 timcs) when the measurement of 
the ancilla give the wrong result and then wc just need to repeat the state preparation. 
This failure does not corrupt the quantum information as the qubtis have not yet inter- 
actcd with these photons. When the measurement of the ancilla give the right vàlues, 
we know that the state is the correct one for teleportation. 

The success probability for this CSign gate is \. Exercise: given that the dashed- 
dotted box in fig. Eümakes a state of the form | (|1010) + |1001) - |0110) + |0101)) 56g7 , 
show that this circuit executes a C-sign between Q\ and Q^- Exercise: what error 
occurs on the information qubits if the detectors R\, R$, i?5 and R-? do not measure the 
correct sequence? 

2'5.7. Incrcasing the probability of success. The task that remains now is to improve 
the probability of success of the teleportation. If we increase the complexity of the state 
produced in the dashed-dotted box of fig. |55l we can increase the teleportation success 
probability asymptotically closc to one. We can generalise eq. I|17f) by increasing the 
number of modes: 

\t n ) = —L=j2\iy\o) n - j \o)\i) n -' 

where the exponent k means that this state is repeated k times, i.e. |1) 2 — > 
Or in the logical basis 
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When we use this preparation state for teleportation the success probability scales as 

n i 1 

n+l — n+l ■ 

2'5.8. Gcncralised beam splitter. When we use the largcr entangled state \t n ) wc 
need to use a generaliscd beam splitter: the n+l point Fourier transform. 

Wc definc the n + l point Fourier transform with F n +\. The matrix elements are 
given by 

■27Tfçi. 

u(F n+1 ) = — 

fci s/n + 1 

Exercise: show that for n = 1 we get back our anti-symmetric 50-50 beam splitter 
transformation. 

As in fig. 1231 we measure all the outputs of the n + l point Fourier transform. If we 
detect k photons and < k < n + 1, then we know the qubit has been teleported to 
mode n + k. If we measure either or n + 1 photons the teleportation has failcd. 

For example, consider the case when n = 2: 



and 




Exercise: Find the linear circuit that implement this unitary transformation. 




Fig. 27. - Teleportation with | Í2 ) . The dashed-dotted box produces |Í2 ) - The probability of 
success of this teleportation is 2/3. 
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Whcn F3 is applied to the incident modes 0, 1 and 2, the modes are transformed as 
follows: 

«0 -> + +4) 

a l ^f( a o + e 3 a l+ e 3 a l) 

„ + 1 /„ + _ÍHÍ»t 27T» . J-v 

«2 ^= K + e 3 «1 + e 3 a 2) 

Now say we want to teleport the state a|01) + /3|10) (where the second mode of this state 
propagates unchanged) with the entangled resource state ^2}- 

Once we apply the 3 point Fourier transform to the states a|01) + /3|10) and | í 2 ) , we 
have 19 mcasurcment outcomcs: Ro,Ri,R 2 = 000, 001, 010, 100, 011, 101, 110, 002, 020, 
200, 012, 021, 201, 102, 120, 210, 300, 030 and 003. We are only interestcd in measuring 
total photon numbers between and n + 1 = 3, that is, only 1 or 2 photons in total. 

The measurement outcomes we are interested in are givcn in the table below, along 
with the mode the information is output in and the probability of that measurement. The 
following identitics will be useful when working through this example: e í 3 7T -\-e~ % 's 7 ' = — 1, 
e if 7T + e -if7r _ _ 1) e i§ ix + 1 = e if an d e -ii* + 1 = e-'f 

The total probability for successful teleportation is 3(|) + 3(-§f) + 3(^y) = |. This 
corresponds to with n=2. As in the case of the previous subsubsection, when the 
teleportation fails we know it and then the state is projected in the Z basis. 



Measurement Output Mode Probability 



|100) ±(a|0) +/3|1)) 3 

|010) \{e^ra\0)+P\l)) 3 

|001) |(e- a f i a|0)+^|l)) 3 

1200) _^( Q | )+/3|1)) 4 

|020) 5^(«|0> +e^m) 4 

|002) ^(a|0> +e- a 5 i /3|l» 4 

|110) _^(_ Q | ) +e ¥ / 3|i)) 4 

|101) _^(„ Q | ) +e -¥^|i)) 4 

1011) ^(-a^-^l)) 4 
Table II. - The telelported states when using | £2 ) as an entanglement resource. 



2'5.9. Bounds on success probabilities. We have shown in this lecture that we can 
perform the NS-i gate operation a|0) + 7|2) — > a|0) + (3\1) — -y[2) with a success 

probability of i (fig. II 1|> . We have also shown that using two of these NS-i gates we 
can perform a CSign with success probability (hg- I12f> . However, it is not known 
whether these probabilities are optimal. 
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To date, the best NS-i gate success probability is as above and the best CSign gate 
succeeds with a probability of ^ ^HJ- The ^ CSign gate is shown in fig. EHl 



7T 




Fig. 28. - A CSign gate with success probability J^. Here 6»i = 54.74°, 4>i = 0, 62 = 
-54.74°, 02 = 0, 3 = 54.74°, <fo = and <9 4 = 17.63°, 4 = 0. 

In this circuit the two qubits Q\ : a|01)i5 + /?| 10) 15 and Q2 : 7|01)26 + 5| 10) 26 are to 
have a CSign applied, where modes 5 and 6 are left unchanged. Exercise: check that 
this circuit implements a CSign on mcasuring single photons in modes 3 and 4. 

Knill has shown that the CSign gate has an uppcr bound of | when using two ancilla 
modes, each with a maximum of 1 photon per mode |17j . Knill has also shown that the 
upper bound for any NS-i gate with one ancilla mode is \ ^Jj. Scheel et. al. conjecture 
that the upper bound for the NS-i gate is j, regardless of the dimensionality of the 
ancilla space |T5] . 

3. LOQC and Quantum Error Correction 

3T. Improving LOQC: beyond state preparation. - In the last section we showcd that 
scalable quantum computing was possible with linear optical elements, single photon 
sources and photo-detection. The problem of making a two qubit gate with probability 
of success arbitrarily close to one was transferred to a state preparation problem. By 
using gcneralised entangled states of the form |í„), the probability of success scaled as 
-^çi- Howcvcr, states of the form \t n ) are complicated to make. To date the best schemes 
need resources that are exponential in n |19l I2Ü| . Instead, we may ask ourselves if we 
can incorporate quantum error correction \2l\ into the LOQC proposal. Is it possible 
to correct for the incorrect measurements in the bàsic linear optical teleportation (fig. 
I2U|> ? Can we use quantum error correction along with smaller |í„) states to increase the 
probability of successful gates giving scalable quantum computing? 
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To answer this we íirst need to know what errors are inflicted on the quantum infor- 
mation when an incorrect measurement is made in the teleportation in fig. 1201 

In the previous section we went through this bàsic teleportation with linear òptics 
and showed that if we want to teleport the state a|01) + /3|10) using the entanglement 
resource state -75OOI) + 1 10) ) , we need to measure 1 and only 1 photon at the photo- 
detection stage (fig. 123)) . The probability of success was \. But what happens when we 
measure or 2 photons? 

Recalling from the last section, our input state is transformed to 

-> - (\/2a|0101) - a|1100) + a|0110) + /3|1001) + /3|0011) - /3|2000) + /3|0020) 

just before the measurement in fig. 1231 (remembering the 2nd mode is untouched). If 
the detectors measure a total of 1 photon {R\ + i? 2 — 1) we can recover the original 
information: 

iïi= 0,^ = 1: Í(of|01>+y9|10» 



2 
1 

If the detectors measure a total of or 2 photons (R\ + R2 = mod 2) we have: 



Ho = = 1 : 5(-oo|01)+ai|10)) 



R 1 =0,R 2 = 0: 4= Q! o|ll) 



R = 2,R 1 =0: — ^=ai|00) 



R 1 =0,R 2 = 2: -^=ai|00) 



V2 

1 

72' 
1 



Note that when we measure or 2 photons we have measured our qubit, collapsing our 
"quantum" information. This means we are free to add a single photon, or vacuum mode. 
When we measure photons, we obtain the state 1 01} = |0) g and when we measure 2 
photons, we obtain 1 10) = When wc measure or 2 photons we effectivcly measure 
our information qubit in the Z-basis. Exercise: check that the probability of mcasuring 
0, 1 or 2 photons in total adds to 1. 

Knowing this we should be able to use quantum error correction to somehow encoded 
our quantum data and then correct when the teleportation fails, i.e. when the measure- 
ment give either or 2 photons in fig. 1231 This means we need to find quantum error 
correcting codes that deal with an error model corresponding with projection in the Z 
basis, where we know which qubit was projected. From this point on wc'll be dealing 
with Z-measurement errors that arc appearing when we try to make a two bit gate on 
the qubits. We will assume that the one bit gates (implemented by linear òptics) are 
error free. But first we need to know what a quantum error correcting code is. 
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3'2. Quantum Error Correcting Codes. - 

3'2.1. What are thcy? . We can encodc our quantum information with codc words 
to protect against errors. Classically, when the error model is given by independent bit 
flips (X), this can be done via the repetition code: — ► 000 and 1 — » 111. 

When we generalise the classical theory to the quantum one, we also have to worry 
about phase errors (Z operators and the combination of bit flip and phase errors: Y 
opcrators). A necessary and sufficicnt condition for a codc with basis code words {IV')} 
to correct for errors {E} is 

(18) (^\ElE b \^) = C ab 8 l3 

where |^), |^-) G {\tp)} and E a ,E b E {E}. 

Say we have errors of the form E\ = i(J + Z), Ei = i(il — Z). What does this do 
to our qubit \ip) = a\0) q + /3|1) ? ? 

EM) = a\0) q 
EM=(3\l) q 

We have measured our qubit in the Z-basis. How can we encode to correct for this error? 
We could try the following encoding: 

(19) |0) = ^=(|00) Ï + |11) Ï ) 

(20) |1> = -^(|01), + |10),) 

If we measure the first mode of our encoded qubit in the Z-basis we have 

-J=( a |00) 9 +a|ll) g +/3|01) g +/3|10) ç ) 

= -^{|0) g (a|0) g + /3|l),)} + -^|l),(a|l) g + J 9|0) g )} 

and we sec that if we measure mode 1 to be a |0) g , we get our original information qubit 
back and if we measure mode 1 to be a |l) q , we can get our original information qubit 
back once we apply a bit flip. Now we need to check that this code satisfies condition II 81 
for the errors E\ and E-^. We assume the errors are on the lst qubit, so E\ = i (l+Z) ® 11 
and E 2 = |(ll- Z) ® 1: 
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<0|£?Í1|0> = i (014110) = | (0|i5 1 t £; 2 |0>=0 (0|^i|0>=0 

(0\E\l\l) = (0|S+1|1)=0 (0|£ÍE 2 |1) = (OI^Sill) = 

(Tj^ljo) = (T|£;|l|0)=0 (T^Sajo) = (Tl^^ijo) = 

(T|SÍ1|Ï> = i (T|4l|T) = | {T|£Í£ 2 |T>=0 (T|StE 1 |ï>=0 

Table III. - Testing condition \18í for the errors E\ and E%. 



But this doesn't teli us how to correct for the errors E\ and E2, only that correction 
is possible. 

3'2.2. Z-measurement Quantum Error Correcting Code (QECC). How can we correct 
for the Z-measurement errors E\ and E<p. There is a formalism called the stabilizer 
formalism (see [32]) that gives us a way to do this. The stabilizer is defincd by a set 
of operators for an n qubit system whose common eigenvectors define a 2 fc -dimensional 
subspace (the code). That is, the stabilizer is defined as the set of operators {Mi} 
that leave the code word space {iV'j)} invariant: Mi\tpj) = \tpj) We will focus on 

operators formed by tensor products of Pauli operators. We have the two code words: 
|0) = ^(|00), + |ll) g ) and |T) = ^(|01), + |10) 3 ). Thcse are stabilized by the two 
operators 11 ® 1 and X ® X. Exercise: check that XX stabilizer the Z-measurement 
code. 

For a stabilizer code we can detect all errors that are either in the stabilizer or anti- 
commute with any member of the stabilizer group. Say we have the QECC with code 
words \i[>j) and Mi is an element of the stabilizer S. And say E is an error such that 
{Mi, E} = 0, thcn MiE^j) = -EM^j) = -E^). So E^j) is an eigenstate of M», so 
measuring Mi will teli us if E has occurred. By measuring all the stabilizer generators 
we idcntify the error syndrome with this mcthod. 

If a QECC can correct for errors E and F, then it can also correct for the error 
aE + bF. From the form of Ei and E2 we need only consider correcting the error Z (g> 1. 
More generally, when we consider correcting for a Z-measurement error on either qubit 
of the code words H19f) , 120|) , we need to consider the errors Z ® II and \ ®Z. Both Z ® H 
and 1 (g> Z anti-commute with X Cg> X. Exercise: check this. 

The operator X g) X has the matrix form: 




/ 

and has eigenvectors: 







X ®X 







( 












1 


(! 


i)« 


9 (i 


j 


= 





1 






V 1 



IIXio,») = (±|00), + |11),)/V2 
\XX llfi i) = (±|01), + |10),)/V2 
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with the following eigcnvalucs: 

X®X\XX 00 ) = -\XX 00 ) 
X®X\XX Q1 ) = -\XXoi) 
X®X\XX W ) =+\xx w ) 
X®X\XX n ) = +\XX n ) 

For example, consider the first qubit of our encoded state = a\0) + j3\l) = 
-±=(a|00) 9 +a|ll) 9 + 0\OÏ) q + (3\10) q ) is measurcd in the Z-basis: 

\^)^±\0) q (a\0) q +P\l) q ) or -L|l) 9 («|l) 9 + /3|0) g ) 
What if we measure the first qubit in the +Z eigenstate |0) g ? becomes: 
i (a\XX w ) - a\XX 00 ) + 0\XX n ) - P\XX i)) 

If we measure X ® X and get +1 we have: ^(a|00} 9 + a|ll) g + 0\Ol) q + /3|10) 9 ), 
our original encoded state. If we measure X ® X and get —1 we have: ^j(a|00) 9 — 
a|ll)q + /3|01) ç — P\10) q ), and we need to perform the Z operation on the first qubit to 
get back our original encoded state. 

Now what if we measure the first qubit in the — Z eigenstate |l) g ? becomes: 

^(a\XX m ) +a\XX 10 ) + 0\XX ol ) + 0\XX n )) 

If we measure X ® X and get +1 we have: ^(a|00) 9 + a|ll), + (3\01) q + f3\10) q ), 
our original encoded state. If we measure X <g> X and get —1 we have: -^=(—a\00) q + 

a\íl) q — /3|01) q + (3\10) q ), and we need to perform the Z operation on the first qubit to 
get back our original encoded state (leaving us with an ovcrall phase factor). This can 
be summariscd in fig. 1291 

Alternatively, a simpler method is shown in fig. I3UI where we have (UV)g = (U <g> 
V)g = cos | — i sin | (UV) and (U ® V) 2 = I. Exercise: show that fig. 1201 is equivalent 
to fig. 1291 correcting for Z-measurement errors on the top qubit. 
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Fig. 29. - Quantum error correction for Z-measurement errors. 




Fig. 30. - An alternate quantum error correction circuit for Z-measurement errors. 
3'3. Properties of the Z -measurement QECC '. - 

3'3.1. State Prcparation. If wc start with the state \%p) = a\0) q + f3\l} q , how can we 
prepare the state \tp) — a\0) + /3|ï)? To do this we need an ancilla |0) g state, two singlc 
qubit rotations and one two qubit rotation, as seen in fig. 1311 

Exercise: what singlc qubit rotations arc needed to makc ZZgg from a ZYgo? Notc 
that depending on the choice of a and /3 we can prepare the encoded eigenstates of X , 
Y and Z. X, Y and Z eigenstate production are error free. Exercise: check that fig. 
I3l·l pcrforms the desired transformation on a\0) q + /3|l) g . 

3'3.2. Singlc qubit rotations. How do we perform encoded X, Y and Z operations? 
It is to verify that the following operators indeed perform the desired operations: 

X = X^ (orI< 2 ') 

y = y(D Z (2) ( or Z (D r (2)) 

z = z^z^ (or -y(Dy( 2 )) 

Exercise: check this. 
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©=1 



IQ 



|o>. 



Fig. 31. - State preparation for the Z-measurement code. 



Any single qubit rotation can be made from Xa, and Zgo (where Ug = cos | — i sin | U 
and U 2 = 1). Exercise: show this by showing how to make Ya, from Xa, and Zgo. 
The operator Xa, is implcmentcd as follows: 



Xa, = cos - - i sin - 



X = cos^-isin^X( 1 )/( 2 )=XÍ 1) 7( 2 ) 
2 2 * 



The operator Zgo is implcmentcd as follows: 

Z 90 = ^=(1 ÏZ) = -^(7 - iZ«Z( 2 )) = (Z^Z^go 

The (Z^Z^^go operation is closely related to the CSign as can be seen in fig. Wc 
will show later that teleportation tcchniques can be used for this gate as it is not error 
frec. 




Fig. 32. - The equivalence between the CSign and (Z {1) Z {2) ) 90 operation. 



The Z±go gates can be implemented from error free (in our present model) single 
qubit rotations. 
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3'3.3. Measurements. After our qubits have been encoded with the code words from 
eqs. (|19[) and (|2t)(l . how do we measure in a particular encoded basis? 

To measure in the encoded Z eigenbasis: measure both qubits in the Z eigenbasis. 
The encoded states |0) and |ï) arc distinguishcd by the parity of the qubits. This can 
be done error free. To measure in the Y basis wc nced to rotate the qubit using a Xqq 
rotation, and then measure in the Z basis as mentioned above. This again can be done 
error free. If we wanted to measure in the X basis, we would nced to first rotate the 
qubits using a Ygo- This rotation. in term of the bare qubits, requires a YZgg gate and 
cannot be done error free, it will require a CSign gate and thus a chance of an error. 

3'3.4. Two qubit rotations. In order to have a universal set of encoded operations we 
need an encoded entangling operation. The (ZZ)go rotation sufhces for this. This can 
be used to make an encoded CNOT, as shown in fig. 1321 



(ZZ) 90 = -j=(l-iZZ) 

= _L(i_ ízV>z®zWzW) 

v2 

= (Z^Z^Z^Z^) 9Q 

The circuit to implement (ZZ)go is given in fig. 1331 Exercise: confirm this. 

Unfortunatcly fig. 1221 does not readily yicld a logical gate with significantly less error 
than using the large entangled states \t n ) shown in the last section. To do this requires 
teleportation techniques, as will be explained in the following subsections. 




Fig. 33. - A (ZZ)go gate between qubits \Q)i and \Q)2- 



3'4. Summary so far. - With linear òptics we can use teleportation along with the 
large entangled states |í„) to perform two qubit operations at the expense of having Z- 
measurement errors. Without correcting for thcse Z-measurcmcnt errors the probability 
seales as — tt. 

n+l 

To correct for thcse errors we use the encoding |0) = "^(lOO)^ + |H) g ) and |1) = 
-^=(|01}q + |10) q ). Encoded X, Y and Z eigenstate preparations are error free as are 
encoded Y and Z measurements. 
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We can perforin encoded operations on this code. X$ is error free however ÍTgo 
requires a two qubit gate so is not. Similarly (ZZ)go is not error free. 

3'5. Threshold for Z -measurement QECC . - 

3'5.1. Accuracy threshold Theorem. An implementation of information processing is 
sealable if the implementation can realise arbitrarily many information units and oper- 
ations with fixed accuracy and with physical resource overheads that are polynomial in 
the number of information units and operations. 

Accuracy Threshold Theorem: If the error per gate is less than a threshold, then it is 
possible to efHcicntly quantum compute arbitrarily accurately. 

This theorem can be proved by designing fault tolerant gates so that the error model of 
the encoded gates is the same as the one of the bare qubits. Then wc need to concatenate 
these gates. If the error per qubit at one level of encoding is be less than the error per 
qubit at the next level of encoding we can keep this concatenation and reduce the effective 
error rate. 

3'5.2. Nice teleportation. In order to understand how to rcach a threshold for quan- 
tum error correction in LOQC we need to first consider the "nice" teleportation circuit, 
shown below in fig. &\ The reason why this circuit is "nice" is because of how it effects 
the quantum information whcn there is an unintentional measurement error. 




Fig. 34. - A nice teleportation circuit. The state produced in the dotted-dashed box is \tX2) = 
i(|00) g +í|01}<, + |10) 9 -i|ll) g ). 

We want to teleport the state (a\0) g + /3\1) g j 1 . The action of the correction gates Cy 
and Cz for each combination of S'í and S2 are given below: 
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Si 5*2 C y C z 

-1-11 1 

+1 -1 Z 11 

-1 +1 11 X 

+1 +1 Z X 



Exercise: confirm that if we measure 



\Y+) 1 \Z+) 2 weget (a\l) q 
\Y+) 1 \Z-) 2 wcgct (a|0), 
\Y-) 1 \Z+) 2 wcgct (a|l), 



/3|0) 9 ) 3 

/?|1) 9 ) 3 
/?|0) 9 ) 3 



|y-)i|Z-) a wcgct (a|0) ? + /3|l) ? ) 3 

bcforc thc corrcction gate. 

If an error occurs during the Zg^ state preparation gate, we just try again as 
they do not corrupt the quantum information. When we consider errors we need to 
analyse the effect of an unintended Z measurement error on qubit 1 and the effect of an 
unintcnded Y measurement error on qubit 2 |23| . 

The motivation behind looking at unintended Z and Y measurement errors can be 
seen by looking at the rclationship between CSign and the ZYqq gate. First consider a 
CSign gate. As we saw in section fe'll when we use teleportation to implement a CSign 
gate it can fail with either a Z-measurement error on qubit 1 or a Z-mcasurcment error 
on qubit 2, as seen in fig. 1351 



OR 



Fig. 35. - Errors that can result from the teleported CSign gate. 



Now consider how these unintentional Z-meansurement errors change when we look 
at the ZYgQ gate. Exercise: using figs. II UI andlÏÏ^lshow that CSign and ZY$q are related 
as shown in fig. Eïïfl 

Tracing back the potential failures of the ZYqq gate due to our error model using 
figure w e get that they can originate from a failure of the CSign gate on the first 
(control) qubit, leading to a projection in the Z basis or they can originate from a failure 
of the CSign gate on the second (target) qubit, leading to a projection in the Y basis. 
And, again, this ineludes thc knowledge of which qubit has failcd (as we know this for thc 
CSign gate). Thus the error model after the teleportation is similar to the error model 
for the unencoded qubits. 
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In more detail, consider the case of a Z mcasurcmcnt error on qubit 1 in fig. 1341 
We need to first work out what the state looks like at point A in fig. 122 

(21) J-í ia \ooi) q + a \010) q + /3|100} 9 - i0\Ul) q ) 

y 2 ^ ' 123 

If we measure qubit 1 in the Z basis we have 

|Z+>i(ía|01>, + «|10) ? ) 23 
or 

|Z-)i(i/3|00>, + /3|ll) g ) 23 

If this error does occur, our quantum information is corrupted and the gate failcd. The 
failurc mode is again a projection in the Z basis. 

Next consider that we have a Y mcasurement error on qubit 2 in fig. 1341 Oncc again 




A 



Fig. 37. - Teleportation with a y-measurement error on qubit 2. 

the state at point A in fig. |ÏÏ7|is given by eq. I|21() . If we measure qubit 2 in the Y basis 
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we have 

\Y+) 2 (ia\0l) q - ia\00) g + p\W) g - /?|H> 9 ) 13 = \Y+) 2 (-ia\0) q + ^|1>«) 1 (|0>« - l 1 )^ 

or 

|y-)a(*a|01>, + ia\00) q + /3\10} q + /9|H),) i3 - |^) 2 (*«|0) ? + /3|1) 9 ) i (|0), + |l) g ) s 

If this error does occur, our quantum information can be recovered and we can re-use 
qubit 1 for teleportation. 

3 '5. 3. Teleportation with error recovery. Now say we have an encoded qubit a|0) + 
that experiences a Z-measurement error on qubit 1. We saw in fig. 1301 that we could 
correct for the error. When we teleport each part of the encoded qubit and then correct 
for the Z-measurement error on qubit 1 we have: 




Fig. 38. - Correcting for Z-measurement error on qubit 1. 

We can commute the (Y^ X^)-go, measurement and (Y^ X^)go gates to the state 
production stage, shown in fig. |2ÜJ 

Exercise: how are the corrections Si, S2, S3 and S4 modified after the {Y^X^ 2 ')-qq 
and (y«X( 2 )) 90 are commuted through? 

We take / to be the probability that teleportation fails. We want to calculate the 
probability of networkl39lto fail: F r . We can assume the state preparation has succeeded 
as it is done offfine. We have to look at the modes of failure for the (Z^Y^)go gate 
(marked with an asterisk, *) after the state preparation. First we have a term for the 
probability that qubit 4 will be projected in the Y basis (which fails with a probability of 
/), which can be undone, and we retry the wholc circuit (retrying fails with a probability 
of F r ) and thus we get JF r . Next we have the probability that the Y projection of qubit 
4 succeeds (1 — /) but qubit 2 is projected into the Z basis (/) and the circuit fails: 
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Fig. 39. - Correcting for Z-measurement error on qubit 1 with the error correcting rotations 
commuted to the state preparation stage. 



(1 — /)/• We therefore have 



(22) 



F r = fF r + (1 - /)/ -» F r = /, 



i.c. the probability of this circuit failing is the same as the probability that the tele- 
portation fails which is the same as the corruption of one of the qubits in the CSign 
gate. 

3'5.4. Encoded Zgo Gate. Remember that the encoded Zgo gate is given by: Zgo = 
(Z (1 )Z (2 )) 90 . As with the error correction after teleportation circuit (figs. 1381 and lïïü|l . 
we can commute the (Z^ Z^)gg back to the state preparation stage, as shown in fig. 

im 

We want to calculate the total probability of failure of recovery (F z ) for this circuit. 
The probability that both teleportations fail is / 2 . That is, the probability that both 
qubits 1 and 2 are projected into the Z basis is / 2 . The probability that one qubit is 
successfully teleported but the other is not, and recovery fails, is (1 — f)f 2 - That is, 
the probability that the top teleportation succeeds and qubit 4 is projected in the Y 
basis and recovery fails with a projection of qubit 2 in the Z basis is (1 — /)/ 2 . The 
probability that one teleportation succeeds but not the other, recovery succeeds and we 
re-attempt, is (1 — f)fF z . That is, the probability that the top teleportation succeeds 
and qubit 4 is projected in the Y basis and recovery succeeds, so we retry the entire 
circuit is (1 — f)fF z . We therefore have 



(23) 



f + (1 - f)f + /(l - f)F a 



F z = 



2/ 2 - / 3 
l-f + f 



and when the teleporation fails, the mode of error is a projection in the Z basis. 
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Fig. 40. - An encoded Z 90 gate: Z go = {Z w Z (2) ) 9 o- 



In order to reduce the probability of failure one would perform the teleportations 
sequentially. To perform (ZZ) 90 = (Z^ Z^ Z^ Z^) 90 we would simply have the samc 
circuit with two more copies of \tX2)- A similar result for F zz is obtained |28| . 

3'5.5. Thrcshold. As mcntioned earlier, the error per qubit for one level of concate- 
nation must be less than that for the next level up. In our case, the error for the ZZgo 
(or Control-sign gate), /, must be equal to the error for the next level of encoding, F zz 
(which is equivalent to F z ). If wc plot F zz vs /, we see that F z = f whcn / = |, as 
can be seen in fig. 1411 Thus to cnsure that the fault tolerant error correction proccdurc 
reduces the error probability we need f < i. 

3'6. Other Errors. - So far we have only considered errors due to incorrect measure- 
ments at the teleportation stage. There are many other types of errors that may also be 
important, such as phase errors, incfficicnt sourccs, mode absorption/dispcrsion, impcr- 
fect beam splitters, ineflïcient detectors, . . . We will consider inefficient detector errors 
in the next subsection. 

3'6.1. Photon Loss. In the original LOQC paper [Ü], a modification to the gate 
teleportation protocol was shown that allowed for the detection of photon loss, as shown 

in fig. na 

Without going through this circuit explicitly, we know where photon loss oceurs 
through the measurement of ancilla modes. An alternate way to consider loss in telepor- 
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Fig. 42. - Teleportation with loss detection. 

tation is to say the qubit bcing tclcportcd is fully "crascd" , undcrgoing the transformation 

m 

E(p) = \( P + XpX + ZpZ + YpY) 

where X, Y and Z are the usual Pauli operators. This can be considered to be the 
dcpolarising channcl conditioncd on perfect information about which qubit is randomised. 
Consider we have two state teleportations with the state preparation modified to perform 
a CSign gate as in fig. 221 

We want to consider what happens if we find photon loss on the top teleporation (for 
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the moment assume there are no Z measurement errors from the teleportation) . The top 
qubit is fully erased, total information is lost. But what about the bottom qubit being 
tclcportcd? 



I'M? { 




Fig. 43. - A CSign gate with one of the correction gates applied at the incorrect time. 



Bccause the top qubit is completely randomised, when we perform the gate corrections 
dependent on the top Bell measurement, the Z gate (in colour in fig. I43J1 will be applied 
at the incorrect time. So any photon loss on the top teleportation translates to a full 
erasure of that qubit being teleported and a Z erasure of the bottom qubit: 

Z{p) = \( P + ZpZ) 

To correct for erasures we can use the 7-qubit code. This has the stabilzer generators: 
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This code can correct for any combination of 1 or 2 erasures and 28 of the possible 35 
combinations of 3 erasures. Some of the combinations of 4 erasures can also be corrected. 
To recover from an erasure of a qubit in an unknown location we would have to measure 
all six generators. 

We always know the location of the erasure, reducing the number of generators to 
measure down to 2. The 7-qubit code can also correct for Z erasures and Z-measurement 
errors. This is because a Z erasure on a particular qubit is equivalent to an unintcntional 
Z measurement on the qubit of unknown outcome: Z + pZ + + Z_,oZ_ = ZpZ), 
where Z± = i(l± Z). In this case only 1 generator needs to be measured. The circuit 
for this recovery is shown in fig. ^] . 
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Fig. 44. - A recovery from a Z erasure or Z measurement error on the top qubit. All other 
qubits must be erasure free in this case. The remaining 3 qubits in the 7 qubit code are not 
shown. 



In order to estimate how good photon detectors need to be, we estimate the thrcshold 
for fault tolerant computation. The gate error thrcshold was found by Silva et. al. |25| 
to be bounded below by a value between 1.78% and 11.5%. 1.78% results from errors 
due to loss and 11.5% results from Z measurement errors in teleportation. 
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4. Conclusion 

In these set of lectures we have seen how it was possible to makc linear optical clements 
with singlo photon sources and detectors bchavc likc a quantum computcr. A fundamen- 
tal idea was to use quantum teleportation to implement gates. If we have probabilistic 
gates, this teleportation can transfer the dimculty of making a gate to the dimculty of 
making a statc. In our casc an extra difficulty was that the teleportation would also 
induce errors, which could be taken care of. The wholc process was streamlined by using 
quantum error correction, the subject of the last section. 

The work on linear òptics raises interesting qüestions about what quantum infor- 
mation is and the origin of its powcr. Linear optical elements by themselves can be 
simulated efhciently on a classical computer, but when we measure and condition further 
gates constructed from linear elements, this is not possible anymore. It is somewhat 
surprising that this condition on classical information becomes so powcrful. 

In these lectures we have not mentioned the experiments that have demonstrated the 
first steps towards implcmcnting LOQC. This is a very active and interesting subject 
that will be left to another set of lectures. 
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